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ABSTRACT 

The inflessence model has recently been proposed in an attempt to explain both early inflation and present day accelerated 
expansion within a single mechanism. The model has been successfully tested against the Hubble diagram of Type la Supernovae, 
the shift parameter, and the acoustic peak parameter. As a further mandatory test, we investigate here structure formation 
in the inflessence model determining the evolution of matter density contrast 8 = Spm/pm in the linear regime. We compare 
the growth factor D(a) = S/a and the growth index f(z) = din S/dlna to these same quantities for the successful concordance 
ACDM model with a particular emphasis on the role of the inflessence parameters (7, zq). We also evaluate the anisotropy 
spectrum of the cosmic microwave background radiation (CMBR) to check whether the inflessence model may be in agreement 
with the observations. We find that, for large values of (7, zq), structure formation proceeds in a similar way to that in the 
ACDM scenario, and it is also possible to nicely fit the CMBR spectrum. 
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1. Introduction 

It is now widely accepted that we live in a spatially 
flat universe undergoing an accelerated expansion and 
made out of ~ 95% dark ingredients which we know lit- 
tle about. On the one hand, observations of the CMBR 
anisotropy spectrum (see, e.g., Ide Bernardis et al~2 000 
|Hanany et al. 2000||Spergel et al. 20031 and |Page 2004| for 
a review) indicates that the total energy density attains 
the critical one so that the universe is spatially flat. On the 
other hand, the SNela Hubble diagram (. Riess et al. 2 004 
lAstier et al. 2006 ) is a clear signature of the cosmic speed - 
up of the universe expansion, hence discarding with a great 
degree of confidence the old standard picture of a matter - 
dominated universe. Finally, the matter power spectrum 
and the clustering properties of galaxies observed in large 
galaxy surveys ( |Pope et al. 2004| ICole et al. 2005(1 point 
towards the existence of dark matter suggesting that its 
density parameter Qm is of the order of 0.3, far lower than 
the SCDM value Qm = F thus stressing the need of a fur- 
ther component to achieve the critical density. When com- 
bined together, this impressive set of observations moti- 
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vates the entrance on the scene of a new player dominating 
the energy budget and driving the accelerated expansion. 
This elusive and mysterious component is referred to as 
dark energy. 

Although the need for dark energy is clear, its nature 
and fundamental properties are completely unknown. The 
simplest candidate is the well-known cosmological con- 
stant A HCarroll et al. 19921 ISahni fc Starobinski 2000jl . 
which perfectly matches a wide range of observations 
flTegmark et al. 20003| |Seljak et al. 2005| >, hence award- 
ing the name of concordance model to the scenario 
based on A and cold dark matter (CDM). Despite 
this impressive success, the ACDM model is plagued 
by serious theoretical shortcomings, thus motivating 
the search for alternative schemes. This has opened 
the way to an overwhelming flood of papers proposing 
different models for explaining the cosmic speed - up 
and the CMBR anisotropy spectrum with propos- 
als ranging from a dynamical A originating from a 
scalar field (dubbed quintessence) rolling down its self- 
interaction potential (see, e.g., IPeebles fc Rathra~2 003 
and fFadmanab han et al. 2003 for comprehensive re- 
views), to unified models of dark matter and dark energy 
such as the Chaplygin gas IjKamenshchik et al. 20011 
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Bili c et ai. 20021 IBento et ai. 2003jl and the Hobbit 
models (Cardo ne et al. 2004(1 . to braneworld in- 
spired scenarios (|Dvali et al. 20001 ILue et al. 2 004) 
and higher order theories of gravity both in the 
metric QCapozziello 2 002 |Capozziello et al. 20031 
Nojiri fe Odintsov 2003| " ICarroll et al. 20041 

and 



^ and the Palatini 

H Vol lick 20031 [ Meng fc Wang 2003| |Flanagan 2004 



lAllemandi et al. 20041 |Capozziello et al. 200*4 

|Amarzguioi et al. 2005| ) formulations. Although radi- 
cally different in their theoretical aspects, all of these 
models are equally viable from the observational point 
of view, thus indicating that better quality data, higher 
redshift probes, or new tests are in order to break some 
of the degeneracies among different models. 

It is worth noting that both current theoretical 
schemes and observational evidences predict that the evo- 
lutionary history of the universe comprises two periods of 
accelerated expansion, namely the inflationary epoch and 
the present day dark energy dominated phase. In both 
cases, the expansion is usually interpreted as the result of 
the presence of a negative pressure fluid dominating the 
energy budget. It is natural to wonder whether a single 
(effective) fluid may indeed be responsable for both peri- 
ods of accelerated expansion. At the same time, this fluid 
should be subdominant during the radiation and mat- 
ter dominated epochs so as not to interfere with baryo- 
genesis and structure formation. While it is quite diffi- 
cult to theoretically formulate the properties of such a 
fluid, it is, on the contrary, clear what its phenomeno- 
logical features are. Inspired by these considerations, 
some of us have recently proposed the inflessence model 
((Cardone et al. 20 05). Based on a suitable Ansatz for the 
dependence of the energy density on the scale factor a, the 
inflessence scenario has been shown to be able to fit the 
SNela Hubble diagram, and also give correct values for 
the shift K ((Bond et al.TMTI IWang fc Mukherejie 2004| ) 
and acoustic peak A ( Eis enstein et al. 2005(1 parameters. 
While this result gives an observational motivation for 
the model, inflessence is also well founded theoretically, 
since it can be interpreted both in terms of scalar field 
quintessence and as an effective model coming from fourth 
order theories of gravity. 

Motivated by these observational and theoretical re- 
sults, we extend the analysis of the inflessence model 
here by investigating structure formation in this scenario. 
Moreover, we also present a preliminary analysis of the 
CMBR anisotropy spectrum. Both these features are stan- 
dard observables in cosmology nowadays, and it is there- 
fore mandatory to check whether the inflessence model is 
able to survive these tests. 

The structure of the paper is as follows. Section 2 
briefly recalls the main features of the inflessence model 
and explains what the roles played by its characteriz- 
ing parameters are. In Sect. 3, the evolution of matter 
density perturbations is studied in the linear regime, as- 
suming that the inflessence fluid does not cluster on the 
subhorizon scale, which is indeed the case of most dark 



energy models. Section 4 is dedicated to a discussion of 
how the growth index depends on the inflessence parame- 
ters and the constraints that could possibly be extracted 
from a precise determination of this quantity. The CMBR 
anisotropy spectrum is evaluated in Sect. 5, while a sum- 
mary of the results and the conclusions are presented in 
Sect. 6. 

2. The inflessence model 

The key ingredient is the following Ansatz for the inf- 
lessence energy density : 



p(a) = Afa 



1 



a 7 v ; 
a 



1 



a Qj 



(1) 



with a normalization constant AT, slope parameters (/?, 7), 
and two scaling values of the scale factor a/ << a,Q. For 
later applications, it is convenient to rewrite Eq.Q in 
terms of the redshift z = 1/a — 1 (having set do = 1 with 
the subscript denoting henceforth quantities evaluated 
at the present day, i.e. z = 0) : 



p{z) = 7V(1 + zf 1 + 



having defined : 
zi = \ jai - 1 , 
Z Q = l / a Q - 1 • 



i + : V7. , 1 + zq^ 1 



1 + Zl 



l + z 



(2) 



(3) 
(4) 



From Eq.Q, it is quite easy to see that the energy den- 
sity of the inflessence fluid scales like that of dust matter 
(p ~ a~ 3 ) in the range a/ << a << oq, so that, given 
typical values for (aj,ag), the fluid follows matter for a 
large part of the universe history, while it scales differently 
only during the very beginning (a << ai) and the present 
(a >> ag) periods. Moreover, choosing f3 = —3, the fluid 
energy density remains constant for a << a/, thus behav- 
ing like the usual cosmological constant A during the early 
epoch of the universe evolution. Finally, the slope param- 
eter 7 determines how the fluid energy density scales with 
a in the present epoch. 

It is still more instructive to look at the equation of 
state (EoS) w = p/p, where p is the fluid pressure. Using 
the continuity equation : 



p + 3H{p + p) =0 



(5) 



with the Hubble parameter H = a/a and inserting Eq.ljT]) 
into Eq. © , after some algebra we get : 



a fi+z 

3 {2 + Z + Z! 



1 

3 V2 



Z Q 



Z Q 



(6) 



It is worth noting that w does not depend either on 7 
or on zq for high values of z, which is to be expected 
looking at Eq.|j|JJ). On the contrary, these two parameters 
play a key role in determining the behavior of the EoS 
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over the redshift range (0, 100), which represents most of 
the history of the universe (in terms of time) . 

The role of the different quantities (/3, 7, zj, zq) is bet- 
ter understood considering the asymptotic limits of the 
EOS. We easily get : 

lim w{z) = - , (7) 

z — >oo o 

which shows that setting (3 = —3, the fluid EoS approaches 
that of the cosmological constant, i.e., w\ = —1, in the 
very early universe. In general, if we impose the constraint 
< — 1, we get a fluid having a negative pressure in the 
far past so that it is able to drive the accelerated expan- 
sion occurring during the inflationary epoch. It is there- 
fore clear that zj controls the transition towards the past 
asymptotic value, in the sense that the larger is z with re- 
spect to zj, the smaller is the difference between w(z) and 
its asymptotic limit /3/3. This consideration suggests that 
zi has to take quite high values (indeed, far greater than 
10 3 ) since, for z >> zj, the universe is in its inflationary 
phase. 

In the asymptotic future (i.e., z — > — 1), we get : 

lim w(z) = — — (8) 
z— >— 1 3 

so that the slope parameter 7 determines the future evolu- 
tion of the universe. For instance, for 7 = 3, the universe 
finally ends in a de Sitter state (as for the concordance 
ACDM model), while a Big Rip occurs for 7 > 3, as in 
phantom models IjCaldwell 20021 ICaldwell et al. 2003jl . 

Let us now consider the present day value of w that 
turns out to be : 



P 7 



f 1 - 






n- 




V2- 


-zq) 




V2- 





(9) 



where, in the second line, we used the fact that zi is very 
large. Given that zq > 0, to have a present day accelerated 
expansion, wq should be lower than —1/3 so that we get 
the constraint 7 > (2+zq)/(1+zq), Moreover, depending 
on the values of 7 and zq, wq could also be smaller than 
wa so that we may recover phantom - like models. The 
parameter zq then regulates the transition to the dark 
energy -like dominated period. 

In summary, the inflessence fluid with energy density 
and EOS given by Eqs.Q and iJBJ is able to drive the 
accelerated expansion of the universe during both the in- 
flationary epoch and the present day period. Therefore, 
such a fluid plays the role of both the inflaton and 
the quintessence scalar field, hence the motivation for 
the name inflessence (contracting the words inflationary 
quintessence) . 

A general comment is in order here. Although origi- 
nally proposed as a phenomenological Ansatz for the en- 
ergy density, the inflessence model could also be seen as an 
effective parametrization of the scale factor during the uni- 
verse's expansion, which translates into the above scaling 



of p{z) given the known behavior of the matter and radi- 
ation components. This equivalent representation is par- 
ticularly interesting during the inflationary epoch. Indeed, 
if the inflessence fluid plays the role of the inflaton field, 
one should wonder whether its self - interaction potential is 
able to give rise to reheating. Discussing this issue is out- 
side our aims, but we stress that considering the model 
as a parametrization of the scale factor rather than the 
energy density makes it possible to escape problems with 
reheating. 

It is worth noting that, since p scales with a as the 
dust matter energy density for a long period of the uni- 
verse history, the coincidence problem is partially allevi- 
ated. Indeed, the dark energy and the matter components 
track each other for a long period so that their near equal- 
ity today turns out to be a consequence of the relatively 
recent change of the scaling of the inflessence energy den- 
sity with a. However, there is still a certain degree of fine 
tuning since the model parameters have to be set in a 
suitable way so that the transition from decelerated to ac- 
celerated expansion takes place only recently. Moreover, 
although wo < —1 is possible depending on the values of 
(7, zq), the possibility to avoid the Big Rip still remains 
if 7 < 3, although such low values seem to be disfavored 
by the fitting to the dimensionless coordinate distance of 
SNela and radiogalaxies. 

As a final important remark, let us stress that, al- 
though phcnomcnologically inspired, the inflessence model 
may also be theoretically well-founded. Indeed, as shown 
in Cardone et al. (2005), one could obtain Eq.Q as a re- 
sult of scalar field quintessence with a self - interacting po- 
tential, which can be well approximated as a simple power 
law with negative slope. As an alternative explanation, 
the inflessence model may also be recovered as the effec- 
tive fluid like description of a fourth order f(R) theory in 
which the Einsteinian gravity Lagrangian f(R) oc R is re- 
placed by a suitably reconstructed f(R). From a different 
point of view, this also offers the possibility of consider- 
ing the inflessence scenario as an analytical parametriza- 
tion for studying a wide class of diverse models. As such, 
investigating the growth of structures for the inflessence 
fluid gives indirect constraints on these models too. To this 
end, it is worth stressing that we are assuming that the 
inflessence model may be recovered as a particular case 
of scalar field quintessence so that the gravity Lagrangian 
is the standard one, and we can therefore resort to the 
usual perturbation theory for our analysis of the growth 
of structures. Should we choose the interpretation of the 
inflessence model in the framework of f(R) theories, we 
should adopt a different formalism. 

3. Linear growth of fluctuations 

The inflessence model has been successfully tested against 
the SNela Hubble diagram, also taking into account the 
shift parameter (which essentially gives the distance to the 
last scattering surface) and the baryonic acoustic peak 
parameter. Although successfully fitting this dataset is 
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mandatory for any realistic dark energy model, such a test 
only probes the dynamics of the background cosmology. 
Indeed, dark energy also has impacts on the evolution of 
density perturbations. 

As is well known, the universe is homogenous and 
isotropic only on the largest scales. As a consequence, 
while one can still use the standard FRW description when 
considering the dynamics of the universe on the scales of 
interest, smaller scale evolution must take into account 
the inhomogeneities of the spacetime. Using the equations 
of motion for this perturbed metric, one can derive the 
growth of density perturbations. Moreover, it is possible 
to demonstrate that, because of its high sound speed, dark 
energy only clusters on scales that are far larger than those 
of galaxies and galaxy clusters. As a consequence, dark en- 
ergy affects the structure formation process only because 
of its background energy density, which concurs to de- 
termine the expansion rate. This is indeed also the case 
for the inflessence model we are considering here, so that 
we may resort to the standard theory to investigate the 
growth rate of matter perturbations in the linear regime. 

Denoting the matter density contrast with 5 = 
Spm/pmi the perturbation equation reads: 



S + 2H5 - 4nGp M $ = . 



(10) 



It is convenient to change the variable from t to the scale 
factor a so that Eq. ijTU|) may finally be rewritten as: 



5" + 



(InE 



2\n 



5' 



3f2 



M 



2E 2 a 5 



5 = 0, 



(11) 



where the prime denotes the derivation with respect to the 
scale factor a and we have used AmGpu = (3/2)f2j\/iJo a_3 
and defined E 2 = H 2 /H 2 . To study the evolution of per- 
turbations in the inflessence scenario, we have only to 
insert into Ea. lll|) the corresponding expression for E 2 , 
which reads : 

E 2 (a) = H 2 {a)/Hl = f^a" 4 + fi M a~ 3 + flxg(a) , (12) 

where fi r , Qm, and fix are the present day values of the 
density parameters for radiation, dust matter, and inf- 
lessence, respectively, and : 



9(a) 



l + aj/ay fl + a/aQ 



1 + a/ 



1 + 1/ac 



(13) 



Note that, consistent with the position of the first peak 
in the CMBR anisotropy spectrum, we have assumed a 
spatially flat universe so that fl r + fl^ + fix = 1, al- 
though some slight deviations from spatial flatness are 
still allowed by the data when a time - varying dark en- 
ergy equation of state is used. 

Since for a matter -only universe 5 oc a, it is useful for 
studying the effect of dark energy to divide this behavior 
out and switch to the growth variable D = S/a. Starting 



from Ea. Hll(l . it is quite easy to determine the equation 
governing the evolution of this latter quantity 1 : 



D" 



(InE 2 



n 



D' 



M 



3 

a I 1 -lE-tr' 



(\yiE 2 )' 



D 

— = 

a 



(14) 



Equation i|14|) may be solved analytically only in very spe- 
cial cases (see, e.g,. IPercival 20051 and references therein), 
while, for our model (and indeed for most of dark en- 
ergy models), we have to resort to numerical integration 
using the boundary conditions D = and D = 1 as 
a —> 0. Actually, it is not necessary to integrate from 
a = 0, but one may set a — a^s = (1 + zls)~ 1 as the 
initial condition, z^s being the redshift of the last scat- 
tering surface that we compute using the approximated 
relation in Hu & Sugyiama (1996). To this end, we set 
u>b = flbh 2 — 0.0214 in accordance with the nucleosyn- 
thesis constraints ( Kirk man et al. 2003J1 and h — 0.664 
consistent with the Hubble diagram of low redshift SNela 
flDaly fc Djorgovki 2004] ). 

Rather than looking at D{a) directly, it is more inter- 
esting to consider the quantity AD = 1 — D/Da, which 
represents the percentage deviation of the growth factor 
for the inflessence model with respect to that for the con- 
cordance ACDM one. Figure 2] shows AD (multiplied by 
100 for sake of clarity) as a function of the scale factor 
a for different combinations of the inflessence parameters 
(7, z Q ), having set (fi M ,fV) = (0.28, 9.89 xl0~ 5 ) (for both 
the inflessence and the ACDM model) and fixed (J3,zi) 
to their fiducial values (—3,3454) (Ca rdone et al. 20 05). 
Note that these latter parameters play a negligible role 
in our analysis since they mainly affect the evolution of 
the fluid in the very early inflationary epoch. It is worth 
stressing that setting zi — 3454 does not at all mean that 
we are assuming that inflation took place for z near this 
value. On the contrary, as could be easily checked, the 
universe undergoes inflation only for z » zi so that the 
exact value of this latter parameter does not set the end 
of any inflationary period, which could lead to possible 
problems with nucleosynthesis. 

Not surprisingly, the evolution of the growth factor 
highly depends on the values of the parameters (^,zq), 
and both negative and positive deviations from the growth 
factor in the ACDM model may be obtained. Nevertheless, 
some general results may be inferred. First, we note that, 
although deviations as large as 20% may be obtained, for 
most of the parameter space (7, zq) the growth factor of 
the inflessence model is comfortably similar to the ACDM 
one over the range 0.5 < a < 1, i.e., < z < 2, where 
structure formation mainly occurs. Although detailed nu- 
merical simulations should be performed, this preliminary 
result makes us confident that the assembly of galaxies 

1 Actually, one may also use In a instead of a as expansion 
variable. See Linder (2005) for different equivalent equations 
to determine D(a). 
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Fig. 1. AD vs. a for models with 7 from 3.5 to 8.5 (in 
steps of 1) from the top left to the bottom right panel, 
and zq from 1.0 to 5.0 (in steps of 1) from the uppermost 
to the lowermost curve. The other parameters are set as 
explained in the text. 



and clusters of galaxies should have taken place in a way 
that is quite similar to the one in the ACDM model. 

FigureHalso shows that, for a fixed a, the behavior of 
AD with zq depends on what the value of 7 is. For in- 
stance, a better agreement with the ACDM model predic- 
tion is achieved for higher zq if 7 < 4.5, while the opposite 
is true for 7 > 7.5. Actually, AD turns out to depend only 
weakly on zq for 7 > 4.5 so that it is this latter parameter 
that mainly determines the behavior of the growth factor 
D with a in such a regime. To better investigate how AD 
depends on (j,zq), it is therefore interesting to look at 
the contours of equal AD in the (7, zq) plane, which are 
plotted in Fig.EJfor some representative values of the red- 
shift z. Consider, for instance, the results for z = 0.15 (top 
right panel in Fig.EJ. To have \AD\ < 5%, larger values of 
zq are markedly preferred only if coupled with low values 
of 7, while zq ~ 3 is allowed, provided that 7 stays in the 
range (3.7,4.7). As a general rule, the lower the value of 
zq is, the 7 higher must be to still have |AD| < 5%. With 
this caveat in mind, we note, however, that, unless one 
chooses zq > 5 (which is rejected by the SNela fit), 7 < 3 
is disfavored by the requirement that the evolution of den- 
sity perturbations in the inflessence model mimicks that 
in the ACDM one within 5% over the range < z < 1. In 
particular, remembering Eq.©, we argue that models in 
which the universe ends with a Big Rip are preferred. 

It is interesting to note that fitting the SNela Hubble 
diagram with priors on the shift and acoustic peak pa- 
rameters gives 3.17 < 7 < 5.86 and zq < 5.3 at the 95% 
confidence level l|Cardone et al. 2005(1 . Although the best 
fit values (j,zq) = (3.73,0.1) are likely to be excluded 




2.75 3 3.25 3.5 3.75 4 4.25 4.5 




Fig. 2. Level contours for AD in the plane (7, zq) at 
z = (top left), 0.15 (top right), 0.35 (bottom left), 
and 1.0 (bottom right). Contours are plotted for AD = 
(±0.5%, ±1%, ±5%) with short - dashed, solid and long- 
dashed lines, respectively. Note that, for z = 1.0, the de- 
viations are so small that the contour line at AD = 5% 
lies outside the plot. 



because of the large values of AD, it is nevertheless possi- 
ble to find values of (7, zq) that make it possible both to 
fit the kinematic data and give rise to an evolution of the 
structure as similar as possible to that in the concordance 
ACDM model. 

It is worth noting that the above results could be qual- 
itatively explained considering the properties of the inf- 
lessence fluid. As explained in the previous section, the 
energy density tracks that of matter (i.e., p oc a~ 3 ) for 
zq « z « zi, while, for z << zq, p cx (1 + z) 3 1 . 
A large value of zq means that the tracking of matter is 
achieved only at high redshift when the dark energy den- 
sity has become negligible with respect to that of matter. 
For this same reason, high values of 7 are preferred since it 
allows px to increase with z at a slower rate with respect 
to Pm- As a general rule, the preferred values of (7, zq) are 
those that render the inflessence energy density negligible 
with respect to the matter density during the structure 
formation epoch. This is the same mechanism achieved in 
the concordance ACDM model, thus explaining the shape 
of the contours of equal AD in the (7, zq) plane. 

The growth factor D could also be normalized to the 
present day value for the ACDM model, i.e., by set- 
ting D(l) = Da(1). Such an approach may be moti- 
vated considering that structure formation is in remark- 
ably good agreement with the observations on the low 
redshift large - scale structure of the universe. By normal- 
izing to this model at present, we may better investigate 
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Fig. 3. Same as Fig.HJ but now D is normalized to D\ at 
the present day and the curves in each panel refer to zq 
from 1.0 to 5.0, from the lowermost to the uppermost. 



Qm, 7, zq 


a 


di 


d 2 


d 3 


0.28,3.7,0.1 


0.1355 


-0.5725 


-0.0471 


-0.0124 


0.24,5.0,4.0 


0.3371 


-0.6510 


-0.0878 


-0.0367 


0.24,5.5,5.0 


0.3474 


-0.6275 


-0.0814 


-0.0377 


0.24,6.0, 7.0 


0.3513 


-0.5983 


-0.0731 


-0.0372 



Table 1. Values of the fitting parameters in Ea. l|15(l for 

some representative sets of (Qm,j,zq). The first row 
refers to the best fit to SNela fit, while the remaining 
rows refer to models providing a good fit to the CMBR 
anisotropy spectrum. 



how the growth factor deviates from the ACDM one in 
the past. Figure El shows AD (as defined above) as a 
function of a for different values of the model parame- 
ters (7, zq). Comparing Figs-^andEl we immediately see 
that now AD is positive (i.e., D(a) < D\) over the full 
range explored, whatever the values adopted for (7, zq) 
are. Moreover, AD is an increasing function of both 7 and 
zq for a given a, although the dependence on zq turns out 
to be quite weak for larger values of 7. To better investi- 
gate how AD depends on the model parameters, we plot 
the contours of equal AD in Fig.0] for some representa- 
tive values of the redshift z. As expected, at low redshift, 
the deviations are so small that a large part of the pa- 
rameter space gives rise to a growth factor consistent (to 
well within ~ 5%) with the one for the ACDM model, 
so that it is likely that structure formation takes place in 
the same way. Not surprisingly, such a region shrinks as 
the redshift increases as AD gets higher. Nevertheless, at 
z = 1000, AD < 15% over almost the full parameter space 
considered. 
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Fig. 4. Level contours for AD in the plane (7, zq) at 
z = 0.15 (top left), 0.35 (top right), 1.0 (bottom left), 
and 1000 (bottom right). Contours are plotted for AD = 
(1%,2%,3%) with short-dashed, solid, and long-dashed 
lines, respectively, in the first three panels, while AD = 
(5%, 10%, 15%) are used for z = 1000. 

Although an analytical solution of Eq. (|14fl is not avail- 
able and numerical integration is straightforward, we be- 
lieve it is useful to have an approximate expression for 
D(a) to be used in data fitting. To this aim, we have in- 
tegrated Eq. I|14f) for a^s < » < 1, and found that a very 
good approximation is given by the fitting formula : 



D{a) ~ doe"" (1 + d lV + d 2 rf + d 3 if) 



(15) 



with rj — In a, do — D(a = 1), and (a, d\, d 2 , ^3) con- 
stant parameters depending on (fijf, 7, zq). Equation l|15fl 
works very well with an rms error smaller than 2% up to 
z ~ 400 and less than 8% at z = 1100. Unfortunately, 
we have not been able to find satisfactory approximated 
formulae for (a, d\, d 2 , d^) in terms of (f2jy, 7, zq), al- 
though we have generated tables (available on request to 
the authors) for 0.05 < M < 0.55, 2.5 < 7 < 7.5 and 
0.1 < zq < 7.1, which may be easily interpolated to get 
the corresponding fitting parameters. As an example, we 
report their values in Table 1 for some interesting cases. 

4. Growth index 

A quantity that can be measured by the galaxy correlation 
function or the peculiar velocities is the so - called growth 
index defined as : 



f = 



din <5 
din a 



a dd 
S da 



(16) 



From a theoretical point of view, to estimate / for a given 
dark energy model, one may solve Eq. ijTTJl to get 8 = 5(a) 
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Fig. 5. Af vs. z for models with 7 = 3.5 (top left), 4.5 
(top right), 5.5 (bottom right), and 6.5 (bottom right) 
and zq from 1.0 to 5.0 (in steps of 1) from the uppermost 
to the lowermost curves. The other parameters are set as 
explained in the text. 
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and then straightforwardly compute /. However, since we 
often have to deal with numerical integration, it is better 
to directly solve an equation for /, to avoid propagating 
the numerical errors. To this end, one should simply use 
definition l|16fl of / and multiply Eq. Ijllfl by a/S to finally 
get the evolution equation for the growth index : 



f 2 

/' + - + 



(InE 



2y 



/ 



M 



(17) 



Not surprisingly, Eq. 117|) must be integrated numerically, 
using f(aLs) = 1 as the initial condition, and Eqs. (|12fl 
and (|13|) for the particular case of the inflessence model. 

The 2dFGRS collaboration measured the position and 
the redshift of over 220000 galaxies and, from the anal- 
ysis of the correlation function, determined the redshift 
distortion parameter f jb with the bias parameter b quan- 
tifying the difference between the galaxies and dark haloes 
distributions. Using the estimated / jb and the two differ- 
ent methods employed by Verde et al. (2001) and Lahav 
et al. (2002) to determine the bias b, one may estimate 
/ = 0.51±0.1 or / = 0.58±0.11 at the survey effective 
depth z = 0.15. 

Both of these estimates are in very good agreement 
with what is predicted by the ACDM model, so that it is 
interesting to compare the behaviour of / predicted by the 
inflessence model with that of the concordance scenario. 
To this end, we define Af = 1 — ///a, which gives the 
percentage difference between the predictions of the two 
models. This is shown in Fig.0 multiplied by 100 for the 
sake of clarity, considering different values of 7 and zq 
and setting the other parameters as in Sect. 3. It is worth 
noting that Af is always negative, i.e. the growth index / 
of the inflessence model is larger than the ACDM one over 
the whole parameter space (7, zq). It is therefore manda- 
tory to directly compare f(z — 0.15) with the observed 
value (which is in good agreement with the concordance 
model predictions) to check whether the overestimate of 



Fig. 6. Level contours for Af at z = 0.15 in the plane 
(7, zq). Contours are plotted for Af(z = 0.15) = —5% 
(short - dashed) , —10% (solid), and —15% (long - dashed) . 
Models lying to the left of the dotted line have wo > — 1, 
while the ones to the right behave today as phantoms, i.e., 
w < -1. 



/ may be troublesome thus allowing us to put constraints 
on the parameter space. 

As in the case of AD, the dependence of Af on the 
model parameters (7, zq) is particularly involved, so that 
it is better to look at the contours of equal Af in the 
plane (^,zq). These are shown in Fig-IHI where we have 
set z = 0.15 to compare with the observed value of / 
measured by the SDSS survey. It is worth noting that 
such a plot could be used to constrain (7, zq) by requir- 
ing that Af(z = 0.15) be lower than a given threshold 
dictated by the estimated /. It is clear that such a method 
has the potential to severely narrow the region of the pa- 
rameter space (7, zq) in agreement with the observations. 
Moreover, a comparison of Fig^with the projected like- 
lihood contours from the SNela Hubble diagram fitting 
shows that they are orthogonal, so that a combined anal- 
ysis may place strong constraints on the inflessence model 
parameters. Unfortunately, this method is still not appli- 
cable at the moment because the estimated / is still af- 
fected by a large percentage error (~ 20%), which makes 
the test useless, since, as Fig.El shows, Af < 15% over a 
wide region of the (7, zq) plane. This is essentially due 
to the low redshift tested (z = 0.15), but extending the 
measurement to higher z (also with the same percentage 
error) could significantly improve the efficiency of such an 
analysis. Nevertheless, it is worth noting that the allowed 
region of the parameter space (7, zq) lies to the left of the 
w = — 1 line so that phantom like models are excluded. 
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Such a result is nicely consistent with the constraints from 
the SNela fit that also point towards this conclusion. 
Using the definitions of D and /, one easily gets : 



f = D[l 



1 dlnD 
D din a 



(18) 



so that / not only depends on how density perturbations 
evolve, but also on their logarithmic rate of evolution. As 
such, / is a more subtle quantity depending on both the 
expansion history and the structure evolution. As a gen- 
eral remark, we note that Af approaches null values more 
quickly than AD as a consequence of its logarithmic na- 
ture. It is nevertheless interesting to compare the contours 
of equal AD with those of equal Af to see whether nar- 
rower constraints on the model parameters (7, zq) may be 
obtained by imposing the same threshold on both quan- 
tities. Indeed, although the contours turn out to be par- 
allel, they are also shifted towards each other so that the 
overlapping region is significantly narrower than those se- 
lected by the constraint on Af alone. Should we have an 
observationally motivated constraint on AD (as that on 
Af) to the 5% level of precision, we could thus efficiently 
constrain the inflessence parameters (7, zq). 

Indeed, it is likely that future measurements (using 
a larger redshift survey observing more galaxies) should 
lessen the error on the observed f(z = 0.15) to the 5% 
level. In such a case, it would be useful to have an ap- 
proximated formula for f(z = 0.15) as a function of the 
inflessence model parameters. This is given as : 



/(* = 0.15) 



with : 



t 0/1+/2 In 7-I-/3 In zq 
J0" M 



(19) 



(/c/i./a./a) = (0.8778,0.4614,-0.1839,-0.0323) , 

which works quite well (with an rms error less than 
1.5%) for 0.2 < n M < 0.4, 3 < 7 < 6 
and 1 < zq < 7. It is worth noting that, for 
the ACDM model, one has the approximated for- 
mula /a(z) ^ [n M (l + z) 3 ]°' 55 (|Silveira fc Waga 1994| 
|Wang fc Steinhardt 19981 ILokas et al. 2004jl . For the in- 
flessence model, we find a similar formula, but the value 
of the exponent is determined by the model parame- 
ters (7, zq). Note, however, that Eq. (|19|) only holds for 
z = 0.15. Although we have not explicitly checked it, it is 
likely that the same formal expression holds over a large 
range in z, provided that the parameters (/o, /1, /2, J^) are 
suitably changed. 

5. CM BR anisotropy spectrum 

Since its discovery by Pcnzias and Wilson (1965), the 
CMBR has played a fundamental role in cosmology. The 
recent precise measurement of its anisotropy spectrum by 
the WMAP collaboration < |Spergel et al. 2003| ) has further 
increased the importance of such an observable in assess- 
ing the viability of any cosmological model. Unfortunately, 
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Fig. 7. The CMBR anisotropy spectrum for the inf- 
lessence model with zq = 0.1 and 7 = 3,4,5,6,7 (from 
bottom to top). Other parameters are set as explained in 
the text. Data points are the WMAP measurements. 

the large number of parameters entering the determina- 
tion of the anisotropy spectrum makes it quite difficult to 
extract constraints on a given model's parameters from a 
time - expensive likelihood analysis (typically based on a 
Monte Carlo Markov Chain exploration of the wide pa- 
rameter space). This is also the case for the inflessence 
model, so that we will only investigate how the spectrum 
changes as a function of the main model parameters, i.e., 
(7, zq). In the following analysis, we therefore set the fol- 
lowing values for the other parameters involved in the 
computation : 



uj b = 0.04 , n 



M 



0.24 , h = 0.70 , t = 0.17 , n = 1 



with the optical depth r and the spectral index (assum- 
ing no running) n. The remaining inflessence parameters 
(/?, zj) are set to their fiducial values (—3, 3454). To com- 
pute the CMBR anisotropy spectrum, we use CMBFAST 
(Selj ak fc Zaldarr iaga 1996), obtaining the results plotted 
in Figs.[3and|Hl where the inflessence parameters (7,2c?) 
has been chosen with the main aim of highlighting the 
dependence on these parameters rather than providing a 
good fit to the data. 

As in the case of the growth factor D, large values 
of both 7 and zq are needed to get satisfactory results. 
Indeed, while the position of the peaks is essentially inde- 
pendent on the parameters (7, zq), their amplitude is an 
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Fig. 8. The CMBR anisotropy spectrum for the inf- 
lessence model with 7 = 3.73 and zq = 0.1,0.5,1,3,5 
(from bottom to top). Other parameters are set as ex- 
plained in the text. Data points are the WMAP measure- 
ments. 

incresing function of these two quantities. Moreover, at the 
low multipoles (I < 10), the spectrum is significantly over- 
estimated for small values of the inflessence parameters. 
It is worth stressing, however, that very good fits may be 
achieved by suitably tuning the parameters (7, zq). Some 
nice examples are shown in Fig.|3J 

It is worth noting that such large values of (7, zq) stay 
at the upper end of the confidence ranges obtained from 
the likelihood analysis performed in Cardone et al. (2005). 
Nevertheless, they are not excluded, so that we are con- 
fident that a combined analysis of the whole parameter 
space (thus also changing the matter density parameter 
Qm that we have held fixed up to now, the spectral in- 
dex n, and the optical depth r) could pinpoint a nar- 
row region in the parameter space giving a satisfactory 
fit to SNela Hubble diagram, growth index, and CMBR 
anisotropy spectrum. 

6. Conclusions 

The inflessence model has been proposed as a possi- 
ble mechanism to explain both the inflationary epoch in 
the early universe and the present day cosmic speed up. 
According to this scenario, a single fluid with the energy 
density given by Eq. £[J is added to radiation and dust 



Fig. 9. The CMBR anisotropy spectrum for the inf- 
lessence model with {j,zq) — (6,7), (5.5,5), and (5,4), 
from top to bottom. Other parameters are set as explained 
in the text. Data points are the WMAP measurements. 



matter, thus working as the inflaton field at very low a 
(i.e., for z » z{) and as dark energy on the scale a ~ 1 
(namely, for z < Zq). Since this model has been shown to 
be able to nicely fit the SNela Hubble diagram, while also 
giving correct values for the shift and acoustic peak pa- 
rameters, it is worth wondering how structure formation 
takes place. 

To this end, we have investigated the evolution of den- 
sity perturbations in the linear regime, comparing both 
the growth factor D(a) and the growth index f(z) to 
these same quantities in the ACDM model. In particular, 
we have concentrated our attention on the two inflessence 
parameters 7 and zq, which determine, respectively, the 
asymptotic value of the eos (and hence the final fate of 
the universe) and the transition from the matter like to 
the quintessence like scaling of the energy density with a. 
Moreover, since 7 and zq also set the present day value of 
the EoS, one can easily understand that they play a lead- 
ing role in determining both D(a) and f(z). As a further 
test, we have also computed the CMBR anisotropy spec- 
trum for fixed values of the other parameters (especially 
the optical depth r and the spectral index n). 

As a general result, we have found that, using large 
values of 7 and zq, it is possible to work out scenarios 
in which structure formation takes place in quite sim- 
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ilar ways in both the inflessence and the ACDM mod- 
els. Moreover, for these same values, the predicted CMBR 
anisotropy spectrum also nicely agrees with the WMAP 
data. Such large values seem to be disfavored by the fit- 
ting to the SNela Hubble diagram, so that some tension 
between these two different probes is present. However, it 
is worth noting that the constraints coming from SNela 
are rather weak so that it is indeed possible that such a 
conflict is not particularly worrisome. 

It is also worth noting that a precise determination of 
the growth index / (at the 5—10% level) at the low redshift 
typical of present day galaxy surveys or a measurement of 
/ at a higher redshift have the potential to severely con- 
strain the parameters (7, zq). Moreover, such constraints 
are orthogonal to those coming from SNela, so that a joint 
analysis could definitively assess the viability of the inf- 
lessence model and pinpoint a narrow range in the param- 
eter space (fljif, 7, zq). One could also include the CMBR 
anisotropy spectrum in a fully comprehensive likelihood 
test. However, such an approach is likely to be affected 
by strong degeneracies among the five inflessence param- 
eters (Qm, zi, zq) and the other CMBR parameters 
such as the optical depth r, the baryon content Uf,, and 
the spectral index n (and its eventual running dn/dhik). 
To probe such a large parameter space, a Monte Carlo 
Markov Chain approach is mandatory and is left for fu- 
ture works. 

Actually, having determined the growth factor D(a), 
for which we have also found an analytical approxima- 
tion, we may further explore the issue of structure forma- 
tion in the inflessence model. To this aim, one could use 
the estimated D(a) to estimate the critical overdensity 
for collapse at the present day, and as a function of time, 
hence determining the mass function through the Press 
& Schechter formalism l|Press fc Schechter 1 974) for the 
spherical collapse of perturbations, or its generalization to 
elliptical collapse worked out by Sheth and Tormen (1999). 
The mass function is the key ingredient to predicting clus- 
ter number counts, which are known to be a powerful test 
of dark energy models (see, e.g.. lHaiman et al~2 001). In 
many observational applications, it is also interesting to 
check whether the collapsed perturbation is virialized or 
not. As has been pointed out in Percival (2005), dark en- 
ergy also plays a role in this process, and hence it is in- 
teresting to investigate how the inflessence model affects 
this important process. Most of these problems will be 
addressed in a forthcoming paper. 
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